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Abstract 

The relations between the Bernoulli and Eulerian polynomials of higher order and the com- 
plete Bell polynomials are found that lead to new identities for the Bernoulli and Eulerian 
polynomials and numbers of higher order. General form of these identities is considered and 
generating function for polynomials satisfying this general identity is found. 

1 Introduction 

The history of the Bernoulli polynomials B n (s) counts more than 250 years, as L. Euler first studied 
them for arbitrary values of the argument. He introduced the Euler polynomials E n (s), the general- 
ization of these polynomials to the so-called Eulerian polynomials H n (s,p) was made by Frobenius 
[5]. The Bell polynomials [4] also appeared to be useful in combinatorial applications. 

In 1920s N. Norlund [7] introduced the Bernoulli i?i m ^(s|d) and Euler polynomials £4"^(s|d) of 
higher order adding parameters d = {d{}, 1 < i < m. Similar extension for the Eulerian polynomials 
Hn"^{s, p\d) was made by L. Carlitz in [3]. 

The Bernoulli and Eulerian polynomials of higher order appear to be useful for the description of 
the restricted partition function defined as a number of integer nonnegative solutions to Diophantine 
equation d ■ x = s (see [9]). The author showed in [10] that the restricted partition function can be 
expressed through the Bernoulli polynomials of higher order only. 

It appears that there exist relations that connect the Bernoulli, Euler and related polynomials 
and their higher order generalizations to the complete Bell polynomials. These relations presented 
in Section 3 give rise to new identities for the polynomials of higher order, and the identity for the 
Eulerian polynomials is discussed in Section 4. In Section 5 we introduce a more general form of 
the abovementioned identities and find the generating functions for the polynomials satisfying these 
relations. The Section 6 is devoted to derivation of new relations for the Bernoulli polynomials of 
higher order. 
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2 Bernoulli and Eulerian polynomials and their generalization to 
higher order 

The Bernoulli B n {s) and Euler E n {s) polynomials are denned through the corresponding generating 
functions 

4„st 00 -in O&st 00 +n 

e l — 1 ^-^ n\ e l + 1 ^-^ n\ 

n=0 n=0 

Frobenius [5] studied the so-called Eulerian polynomials H n (s,p) satisfying the generating function 

1 00 j-n 

e st ^ = Y, H n(s,p)- (2) 

e l — p n 

K n=0 

which reduces to definition of the Euler polynomials at fixed value of the parameter E n (s) = 
H n (s, —1). The Apostol-Bernoulli [lj and Apostol-Euler [6] are the particular cases of the Eulerian 
polynomials. 

Introduce the generalized Eulerian polynomials through the generating function 



st l-ap ^ u , ^t n 



Y]H n (s,p,a)— . (3) 



e* — p ^— ' ' n! 

n=0 

It is easy to see that 

B n (s) = nH n -i(s, 1,0), # n (s) = H n (s, -1,1), H n (s,p) = H n (s,p, 1). 
Introduce the generalized Eulerian polynomials of higher order H n m \s, p,a\d) as follows 

m 00 



Tl^^ = Y, H tKs,pM^. (4) 



e 1 — Pi i — ' n! 

8=1 rl n=0 



The generating function for the Bernoulli polynomials of higher order B n m \s |d) reads [2]: 



^n^-x^wS- <5) 

i=l n=0 



It can be checked by comparison that 

(n — m) 

the Euler E n m) 

[3] can be written as 



B n m \s\d) = -^-Htl(s,l,0\d), 7n = XI^ 



1=1 



It is easy to see that the Euler En(s\d) [2] and Eulerian polynomials of higher order H n m \s, p\d) 



E^{s\d) = Hl m \s, -l,l|d), (6) 
H(™\s,p\d) = H^(s,p,l\d). 
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3 Bernoulli and Eulerian polynomials of higher order as Bell poly- 
nomials 

The complete Bell polynomials B„(oi, 02, . . .) [5] are defined through the generating function 

expff gi^l^i f. (7) 

\i=l ' / i=0 

Consider Taylor series 

dt ^ ^^tfBiV 



In 



e rfi_l S( ^ / /I" 



where l?j is the Bernoulli number. Using the above expansion in (|5|) we have 

e-exp m-l) 1 " 1 ^! =^B,( S + a 1 ,a 2 ,...)|, (8) 
\ t=i ? / i=o z ' 

where aj = (— l) t-1 Bj<7j/i and er^ = Ofc(d) = denotes a power sum, so that 

#( m )(s|d) =B n (a + oi,02,...). (9) 
Setting here m = 1, d = 1 we find 

_B n (s) = B n (s + ai, a 2 , . . .), «i = — • 

Consider Taylor series for the Eulerian polynomials 

1-p dt 



In 



r 



e M — p 1 — P ~^ 1 — p *! 

where Hi{p) = Hi(0,p) is the Eulerian number. In a particular case with /?j = p in using the 
above expansion we have 

exp ( st ~ r~) exp ( ~ £ ^~^ )<7t | J =X)Bi(s-ffi + oi,a2,...)|, (10) 

^ ^' V i=2 * / i=0 

where a% = —pHi-i(p)ai/(l — p). Comparing (fTOj) to ((7|) we obtain 

i/( m )(s,p|d) =B n (s-a 1 +a 1 ,a 2 ,...). (11) 

The relations (|9llip are the particular cases of a general identity 

Hi m \s,p,a\d) =B n (a- ffi + 01,02,...) (12) 

with Oi = —pHi_i(p, ct)(Ti/(l — ap) and Hi(p, a) = Hi(0, p, a). Note that the value of the polynomial 
is equal to the complete Bell polynomial of the same order with parameters proportional to the value 
of the corresponding polynomial at zero. 
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4 Generalized Eulerian polynomials identities 

Consider a normalized distribution P(t) that can be characterized by the sequence of moments \i n 
or cumulants K n . The moments generation function is P{t) itself and cumulants generation function 
is lnP(t), i.e., 

oo oo 

P(t) = 1 + J2^, lnP(t) = J2^t n - (13) 
n! ^-^ n\ 

n=l n=l 

It is worth to note that in the theory of symmetric functions \x n play role of the complete homoge- 
neous symmetric polynomials, while K n stand for the symmetric power sums. The above definitions 
imply the relation between the moments and cumulants 

n-l 

Pn = K n + V] ( 7 i ) KkfJ-n-k- ( 14 ) 



fc=l ^ ~~ ' 



On the other hand it follows from (fT3l) and ([7]) that 

= B n (Ki,K 2 , . . .), (15) 

so that 

B n («l, « 2) • • •) = K n + ( _ J KfcB n _fe(Kl, « 2) • • •)• ( 16 ) 

fc=l ^ ' 

Using the relation (|12p and the definition ([7|) one finds: 

^ H n m) (s,p,a\d) p ^ fl^(0^Q)ffj , 

ln Z. ^ * =( s -^ t -T^r P ^ ii (17) 

n=0 F t=l 

The last relation leads to the following identity: 

HW (s, P, a\d) = (s- a|d) - E (fc I J) * a ) ff *^-*( s ' a ' d )' 

that is rewritten as 

H™(s,p,a\d) = (s-a^H^^p^ld) - £ ("Wo, p, a)a k+1 H^ k (s, p, a|d). (18) 

Setting here m = 1, di = 1 we find for the generalized Eulerian polynomials 

F n+ i(s,p,a) = (s- l)il n (s,p,a) - — — ^ ( ) #fc(°> p, a)H n _ k (s, p, a). (19) 



5 General case 



Consider a more general problem of construction of polynomials P n (s,r) that can be represented 
by the Bell complete polynomials with parameters aj being proportional to the values of the same 
polynomial at some s = r, that satisfy the generalized version of (I17p 

ln E ^r 1 ^ = <M*. *•)* + <M*. »•) E ^r 1 ^ **» ( 20 ) 

n=0 ' i=l 
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where Qi(s,r) and Q 2 (s,r) are polynomials in s. Then we find 

P n+ i(s,r) = Q 1 (s,r)P n (s,r) + Q 2 (s,r)^2 ( , ) p k{r, r)P n _ fc (s, r). 



(21) 



We want to find the condition on the generating function G(s,r,t) for the polynomials P n (s,r). 
Multiplying the equation (I21j) by t n jn\ and summing over n we obtain 



£ 

n=0 



P n+ i(s,r) 



n=0 



(r, r)P n _fc(s 



n! 



The r.h.s. of the above equation evaluates to Q±(s, r)G(s, r, t) + Q 2 (s, r)G(r, r, t)G(s, r, t), while the 
l.h.s. is obtained by differentiation of G(s,r,t) w.r.t. t, so that 



dG(s,r, t) 
dt 



Ql(s, r)G(s, r, t) + Q 2 (s, r)G(r, r, t)G(s, r, t). 



Consider a particular case of (|22p in the form 

dG(s,r,t] 



dt 



Qi(s)G(s, r, t) + Q 2 (r)G(r, r, t)G(s, r, t). 



The polynomials P n (s,r) generated by G(s,r,t) satisfy the identity 

Pn+l(s,r) = Qi(s)P n (s,r) + Q 2 (r)^( n i )p k (r,r)P n „ k ( S , 



(22) 



(23) 



(24) 



We look for solution in the form G(s,r,t) = H(s,t)F(r,t) that leads to the equation 



dF(r, t) 
dt 



Qi{s) 



dlnH(s,t) 
dt 



F(r,t)+Q 2 (r)H(r,t)F 2 (r,t). 



(25) 



The term in the square brackets in the last equation should be independent of s, so that the following 
condition holds 

Qi(s) _9ta|(M) =M , (t) 

for some function M(t) defined up to arbitrary constant. The condition leads to relation 

H(s,t) = e Ql W- M ® , (26) 
and the equation ()25[) reduces to the Bernoulli equation 

dF(r, t) 



dt 



which has the solution 



M'(t)F(r, t) + Q 2 (r)e Ql(r) *- M(t) F 2 (r, t), 
Qi(r)e M « 



F(r,t) 



The solution of (|23j) reads 

G(s,r,t) = H(s,t)F(r,t) 



CQi(r) -Q 2 (r)e0i«*' 

Qi(r)e Ql ^* 



CQ^r) - Q 2 {r)e^) v 



(27) 
(28) 

(29) 
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where the integration constant C is determined from the condition 

O(a ,r, 0)=Ga (s,r ) = CQ ^ Q2(ry (30) 

Finally we obtain 

G(s r t) = G ( S ,r) Ql (r)e^ 
1 ' ' j Qi(r) + G (s,r)Q 2 (r)(l-eQi(r)t) • ^ 

Setting in the last equation r = 0, Qi(s) = s— 1, Q2(0) = — p/ (1— ap) and Go(s, 0) = (1— ap)/(l — p) 
we obtain the generating function in the l.h.s. of ([3]). 

It is not difficult to generalize the above generating function to the case of polynomials of higher 
order. We define the generating function G(s, r, t, d) for the polynomials of higher order (s, r |d) 
as follows 

m Qi(0)dit 

G(s, r, t, d) = e toM-<W°)]* TT ! — — . (32) 

The polynomial identity reads: 

P#>(.,r|d) = (Qi(.)-Qi(0)(l-<^ (33) 

The case r = s leads to the general solution 

G(s t) = Go(«)Qi(«)« qi( ' )f r34) 

It should be noted that the solutions (|31l34p are valid only when Qi(s) is identically nonzero. 
Otherwise the equation ([22]) reduces to 

dG{S d ^ t] = Q 2 ( s , r)G(r, r, t)G(s, r, t), (35) 
that can be written in the form 

at Q 2 [s,r) 

From the definition of F(r,t) we find G(s,r,t) = exp[Q 2 (s,r)F(r,t)] and obtain G(r,r,t) = 
exp[Q 2 (r, r)F(r,t)]. Thus we arrive at the following equation for F(r,t) 

^^=exp(Q 2 (r,r)F(r,t)), 

which has the solution 

F(M)= > fe <7 )(( r' ) . 

Q 2 {r,r) 

and we obtain the solution of equation (j35|) 

Q2(s,r 



G(s, r, t) = exp 



Q 2 (r,r, 



■hxQ 2 (r,r)(C-t) 



[Q 2 (r, r)(C - t)]-a»(«>»0/O2(r,r) _ (36) 



6 



Setting in ()36[) t = we obtain the integration constant 

J; r>-Q2(r,r)/Q 2 (s,r) 

Q2(r,r 

Finally, we arrive at the solution 



C = —±— G - Q2{r > r)/Q2(s > r) , Go = G(s, r, 0). 



Qi{r,r)/Qi{s,r) 



i(s,r)/Q2(r,r) 



G(s,r,t) = G [1 - Q 2 (r,r)tG 
For r = s equation ([35]) reduces to 

^=%(»)G 2 (M). 
The solution in this case is found from (|34p in the limit Qi(s) — > 

6 Identities for Bernoulli polynomials and numbers 

Using the relation ([9]) and the definition ([7j) one finds: 

inV ^ n g|d V = st + y (-lr 1 ^ f . 

^— ' n! ^— ' i! i 

ra=0 i=l 

The last relation leads to the following identity: 

B^( S \d) = sB^\(s\d) -^( n \-l) k B k a k Btl(s\d). 

k=i ^ ' 

Setting here m = 1, di = 1 we find for the Bernoulli polynomials 

1 n / \ 

B n {s) = sB n ^{s) - - J2 I , I (-l) fe B fc S n _ fc ( S ). 
n fc=i ^ ' 



Finally setting in the above relation s = we arrive at the identity for the Bernoulli numbe 

^n—k — TlB n . 

As all odd Bernoulli numbers (except B±) vanish, one can find 



±{ n \-l) k B k B r , 
k=i v J 



E (f) [1 - (-l) k ]B k a k B^ k (s\d) = -ncnB^Md), 
k=i ^ ' 



from what it follows that 



B^\s\d) = (s-a^Bt^-^^B^Btiisld), 

k=i ^ ' 

1 n / \ 

B n (a) = (a-l)fl n _ 1 (a)--E( jb )^n-*(a) J 

n k=l ^ ' 

-Bn = —B n -\ — — E (j/JB k B n - k , 

71 k=l ^ ' 



where the relation (|47p is the well-known Euler identity and (|46|) is a particular case of the sum 
identity for the Bernoulli polynomials 

S (fy B k(s)B n -k(t) =n(s + t- l)S„_i(s + t) - (n-l)fl n (a + t). 

Rewrite the identity (|45p using the definition of the power sum as follows: 

4 m) (*|d) = (--aijBaWdJ-i^EUJ^SHd) 

1=1 fc=l ^ ' 



k=0 



4 m) (-|d)+Eu P^i-l(^ld) 



8=1 



dud- 



(48) 



i=i 

= (--oiJflJSMd) + ^^( S |d) - I^4-+i) (s | di ), 

where 

di 

The relation (|48[) produces two new identities 

m 

53 4 m+1) (*|d0 = »(« - (rOflJZlCald) - (n - m)4 m) (*|d), 
i=l 

and 

m 

]Tl4 m+1 )( S |d,) = m( s - aOS^L^ald). 

i=l 

Recalling the identity for the Bernoulli polynomials of higher order [7J : 

4 m+1 )(, + di\di) - B^Wdi) = B^\(s\d), 

we have 

m 

5]4 m+1) (^ + di\di) = nsBt\(s\d) - (n - m)B^{s\d). 

i=l 

Consider a particular case of all di = 1, i.e., d = 1; noting that o\ = m we find from ()49p the 
identity derived by Norlund in [7J 



(49) 



(50) 



(51) 



mB^ +l \s\l) = {m- n)B^\s\l) + n(a - m)B i r ^\(s\l). 
It is easy to show that for arbitrary nonzero A the following relation is valid: 

m^ m+1 )(s|A • 1) = (m - n)B^ m \s\X ■ 1) + n(s - Xm)B ( ^\(s\X • 1) 
From the last relation for m = n it follows that 

m 

Bt +1) (s\\ ■!) = (-- Amj^t^lA • 1) = " Ai). 



(52) 



(53) 



(54) 



i=l 
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